Resonant structure of low-energy H3 dissociative recombination 
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New high-resolution dissociative recombination rate coefficients of rotationally cool and hot 
in the vibrational ground state have been measured with a 22-pole trap setup and a Penning ion 
source, respectively, at the ion storage ring TSR. The experimental results are compared with 
theoretical calculations to explore the dependence of the rate coefficient on ion temperature and 
to study the contributions of different symmetries to probe the rich predicted resonance spectrum. 
The break-up energy was investigated by fragment imaging to derive internal temperatures of the 
stored parent ions under differing experimental conditions. A systematic experimental assessment 
of heating effects is performed which, together with a survey of other recent storage-ring data, 
suggests that the present rotationally cool rate-coefficient measurement was performed at SSOtllo K 
and that this is the lowest rotational temperature so far realized in storage-ring rate-coefficient 
measurements on Hj. This partially supports the theoretical suggestion that higher temperatures 
than assumed in earlier experiments are the main cause for the large gap between the experimental 
and theoretical rate coefficients. For the rotationally hot rate-coefficient measurement a temperature 
of below 3250 K is derived. From these higher-temperature results it is found that increasing the 
rotational ion temperature in the calculations cannot fully close the gap between the theoretical and 
experimental rate coefficients. 



PACS numbers: 



I. INTRODUCTION 

The rate coefficient of the dissociative recombination 
(DR) process in Hg has been subject to debate for 
decades [2 [Ij- Over the past few years, there has 
been a convergence in the DR rate coefficient as deter- 
mined by theoretical calculations and measured in differ- 
ent storage-ring and afterglow experiments. The source 
for the remaining descrepancy between the reported the- 
oretical and experimental rate coefficients has been opted 
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to be the internal-state distribution of the Hg ions. It 
has been stated that the H;^ rate coefficient is now suffi- 
ciently well-known for astrophysical models treating the 
H;^ abundance [T]. However, on other astrophysical top- 
ics, such as the reigning H;^ ortho-para spin tempera- 
tures [3] and the observation of the (3,3) metastable state 
d], rotational state-specific rate coefficients are needed. 
Theory and experiment still disagree when it comes to 
collision-energy and rotational-temperature dependence. 
Very little agreement exists in the resonance structure of 
the Hj^ DR rate coefficient, i.e., in the occurence, am- 
plitudes and positions of resonances. Assignment of the 
individual resonances in the rate coefficient of this ele- 
mentary system has until now remained far out of reach. 
Published storage-ring measurements of the rate coeffi- 
cient for rovibrationally cold H;^ report that only the 
lowest rotational states of the vibrational ground state 
are believed to be populated [SHZj. Predictions of cold 
rovibrational rate coefficients differ locally by up to 
an order of magnitude [8 . A recent theoretical investi- 
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gation [S] pointed out that better agreement is achieved 
with the experimental H;^ data if a rotational ion tem- 
perature of 1000 K is assumed. The reason for using a 
higher temperature in Ref. [8^ was to reflect the pos- 
sibility that higher rotational ion temperatures might 
have been present than had been assumed appropriate 
for the experimental conditions. The cold experimental 
H;^ beams could possibly undergo heating during extrac- 
tion, acceleration, and/or circulation in the storage ring. 
No evidence for large heating effects inside the ring that 
could explain a 1000 K or higher temperature has been 
observed to date [H]. Recently, it has been suggested 
that, in the case of ion production with the supersonic 
ion source, the ion-source extraction could have led to 
heating 10] . Further investigations are necessary to shed 
light on the discrepancies remaining between theory and 
experiment in the DR of H;^ . 

In fact, in the long history of DR studies, it has not 
yet been possible to achieve complete quantitative agree- 
ment between theory and experiment for any molecular 
target ion. The closest to achieving full agreement is 
probably the diatomic hydrogen ion, but even that sys- 
tem has not yet exhibited the close convergence of DR 
theory and experiment that has been demonstrated for 
other observables such as photoabsorption. The emer- 
gence during the past decade of a promising theoretical 
candidate for both the mechanism jllj and the quanti- 
tative description of DR [SI [T^HTl] has apparently 
resolved the long-standing uncertainty about the overall 
low energy recombination rate, but the limitations of the 
theory remain unclear. By making a detailed resonance- 
by-resonance test, starting with the present article, this 
will hopefully elucidate both the strengths and limita- 
tions of the theory. It may also prove useful in suggest- 
ing systematic issues that could be affecting experimental 
measurements, such as the ionic rovibronic state distri- 
bution and external field effects. 

We have measured high-resolution DR rate coefflcients 
for rotationally cool and hot H;^ in their vibrational 
ground state for collision energies between and 30 eV. 
The rotationally cool beam was produced with a 
buffer-gas cooled 22-pole trap, while the rotationally hot 
beam was produced with a Penning ion source. The 
use of a cold photocathode electron gun provided us with 
the high resolution |15| . and the greatly improved trans- 
mission of the low-current H;^ beam from the 22-pole trap 
to the storage ring enabled us to gather more statistics 
over Ref.|7]. The rotational ion-beam temperatures are 
investigated through imaging measurements under the 
conditions of the rate-coeficient measurements. We also 
present a re-evaluation of a previous measurement with 
the 22-pole trap [7], which gives an estimate of the ion 
beam temperature produced by the 22-pole trap, shed- 
ding light on a possible heating mechanism. The depen- 
dence of the rate coefficient on ion temperature is inves- 
tigated through theoretical predictions, and the theory- 
experiment comparison permits a first attempt to assign 
the H;^ DR resonance structure. 



II. EXPERIMENT 



The experimental studies were carried out at the stor- 
age ring TSR at the Max-Planck-Institut fiir Kernphysik 
in Heidelberg [161 IE] ■ For the rotationally cool mea- 
surements, the ions were produced in a radiofrequency 
(rf) storage ion source and buffer-gas cooled for 1 ms 
in a cryogenic 22-pole ion trap at a helium density of 
10^^ cm~'^ and a temperature of ~15K. It has been 
shown through spectroscopy measurements and model- 
ing [3 dH] that the cooling of the in the 22-pole trap 
is efficient and produces ions with rotational and transla- 
tional temperatures that agree well with the nominal trap 
temperature. At 15 K, the lowest two rotational states, 
{N+,K+) = (1,1) and (1,0), separated by A = 22.8 
cm~^, of the vibrational ground state are populated, with 
iV+ the total angular momentum of the ion and the 
body-frame projection of the angular momentum of the 
ion. The (1,1) and (1,0) states belong to para and or- 
tho nuclear-spin symmetries, respectively. A differential 
pumping stage between the rf ion source and the rf ion 
trap reduces the partial H2 pressure in the trap chamber 
to 10^^ mbar and minimizes any conversions between the 
ortho and para symmetries. The two states are therefore 
believed to be equally populated instead of Boltzmann 
distributed. For each injection into the ring, up to 10^ 
H;^ ions were stored in the trap, buffer-gas cooled, and 
extracted. The rovibrationally hot Hj]' beam was pro- 
duced in a Penning ion source. 

After extraction, the H;^ beam is accelerated to 4.113± 
0.002 MeV and injected into the ring. Transmission from 
the ion trap to the storage ring has been increased to 
above 50% (as opposed to the < 10% transmission in 
Ref. [7] ) . The transmission of the Penning beam is less 
critical since high currents (fiA) can be produced. When 
present, the vibrational excitation has time to decay in- 
side the ring, leaving rotational excitation only after 2 
s storage time [19 . In two straight sections of the ring, 
named electron cooler and electron target, cold electron 
beams were merged with the H;^ ion beam. The func- 
tion of the electron cooler is to phase-space cool the 
ion beam. Here, the electron and ion beams are veloc- 
ity matched to give 0-eV collision energy at an energy 
spread of 12 meV and an electron density of ^ 10^ cm~^ 
for the entire storage time of 8.5 s. The electron tar- 
get can also phase-space cool the ion beam when set to 
cooling (near 0-eV) energy, but with less efficiency since 
it has a lower electron density of 10^ cm~'^. Primarily, 
the electron target functions as the reaction region to 
study the DR reaction. It provides an ultracold electron 
beam, produced by a nitrogen-cooled photocathode |15| , 
with transverse and longitudinal temperatures of fcTj^ w 
1 meV and fcTy « 30 ^eV, respectively. Two types of 
measurements are described below: rate-coefficient and 
imaging measurements. Because also a re-evaluation of 
and comparisons to earlier measurements are carried out 
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TABLE I: The imaging measurements at near 0-eV energy that are presented and discussed in this paper. The rotational 
temperatures are derived from a fit between simulation and measurement with uncertainties given at a 90% confidence 
interval. The statistical and systematic errors can add up to an uncertainty of about 135 K in the rotational ion temperature. 
The measurements are grouped for each reference. The labels IR refer to the imaging experiments conducted under conditions 
similar to the rate measurements. 



label [Ref] 


ion source 


Troti 




source settings 






cooler 


target 






(AT)stat 


He 




trapping time 


parent 


kTx 


time at eV 


time at eV 






(K) 


(cm-3) 


(cm~3) 


(ms) 


gas 


(meV) 


(s) 


(s) 


IR22p 


22-pole trap 


380 ± 35 


lO^*^ 


10^ 


1 


P-Ha 


12 


0-8.5 


0-8.5 


I22p-09 X 


22-pole trap 


200 ± 75 


10^^ 


10^" 


100 


P-Ha 


13.5 


0-0.5 


0-7.5 


I22n-09 [7] 


22-pole trap 


150 ± 50 


10^^ 


10^° 


100 


n-Ha 


13.5 


0-0.5 


0-7.5 


Ipenn-09 


Penning 


3250 ± 35 








n-Ha 


13.5 


0-0.5 


0-7.5 


Isonic-10 jlOj 


supersonic 


950 ± 100 








P-Ha 


12.5 


off 


0-10 


IRsonic-10 [TD] 


supersonic 


450 ± 100 








P-Ha 


12.5 


0-10 


0-10 



TABLE IL The rate-coefficient measurements presented and 
discussed in this paper. For all rate measurements, the elec- 
tron cooler was set to cooling (near 0-eV) energy for the entire 
storage time. The measurements are grouped for each refer- 
ence. The 22-pole trap pressures of Ref. [6] are similar to 
Ref. [7]- Ref [6] used trapping times of 1-100 ms and Ref. [J 
of 1 ms only. 



label [Ref] 


ion 


parent 


cooler 


target 




source 


gas 


kT_L 


at eV 


at eV 








(meV) 


(s) 


(s) 


R22n 


22-pole trap 


n-Ha 


12 


0-8.5 


0-2 


Rpenn 


Penning 


n-Ha 


12 


0-8.5 


0-2 


R22n-09 [7] 


22-pole trap 


n-Ha 


13.5 


0-11.5 


0-2 


Rpenn-09 


Penning 


n-Ha 


13.5 


0-11.5 


0-2 


R22n-05 W 


22-pole trap 


n-Ha 


13.5 


0-11.5 


0-2 


Rsonic-10 10, 


supersonic 


n-Ha 


12.5 


0-10 


0-2 



here, each performed under different experimental condi- 
tions, an overview of the imaging measurements and the 
rate measurements is given in Tables |l] and |llj respec- 
tively. The respective measurements will subsequently 
be referred to using the labels shown in these tables. 



The internal temperature of the rotationally cool H;^ 
beam has been investigated by analyzing the dissociation 
dynamics of H;^ breaking up into three H atoms through 
an imaging experiment at 0(±0.5)-meV collision energy 
(see IR22p in Table |l]). This measurement was performed 
with both the electron cooler and the electron target at 
near 0-eV energy for the entire storage time. Events with 
three H product atoms were recorded with a position sen- 
sitive detection system as decribed in Ref. [2^. No new 
imaging measurements were performed with the highly 
excited beam. Previous imaging data from the Penning 
source, Ipenn-09, were used instead to estimate an ion 
temperature for the rotationally hot ion beam. 

An improved forward simulation of the dissociations 



was utilized to derive the effective rotational tempera- 
tures from the imaging data. The improvements con- 
tained the following main modifications. The anisotropy 
of the fragment momentum-distribution is now accounted 
for using the Dalitz plot describing the reaction dynam- 
ics presented in Ref. !7j . Additionally, the kinetic energy 
release has been recalculated using the late st findings on 

(see Sec. 



IV I. And finally. 



the ionization energy of Hj^ 
an effective target length approximating the effect of the 
merging and de-merging of the electron beam on top of 
the straight overlap between the ion and electron beams 
has been included. An important difference between the 
current IR22p measurement and the measurements from 
Ref. [7] is that the electron cooler was at cooling en- 
ergy for, respectively, the entire storage time and 0.5 s 
only. To investigate the effect of the electron cooler, the 
measurements from Ref. [7] are re-evaluated using the 
improved simulation. 

The rate measurement, R22n, was performed with the 
electron cooler at cooling conditions for the entire stor- 
age time, similar to the imaging measurement, IR22p. 
This is typical for all rate measurements (see Table |ll]) 
as the electron beam of the target is not permanently set 
to cooling energy as is the case for a 0-eV imaging mea- 
surement. The electron beam of the target was set to 
cooling for the first 2 s only, after which the electron ve- 
locity was detuned to investigate the DR rate as function 
of coUison (detuning) energy. The current rate measure- 
ments cover the range Ed = ~ 30 eV. The count rate 
of the neutral fragments produced by the DR reactions 
in the target was measured with a large surface barrier 
detector (SBD) with an area of 10x10 cm^ situated at 
~ 12 m downstream from the center of the electron tar- 
get. Simultaneously, the charged products as well as 
the neutral H or II2 products from collisions of H;^ with 
residual gas in the ring were recorded with a scintilla- 
tor and the SBD, respectively, as a measure of the ion 
beam intensity. The general detection and normalisation 
procedure is described in detail in Ref. [10]. Here, the 
reference signal was recorded at a collision energy of 10.8 
eV. No absolute measurements were performed. 
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III. THEORY 

The general theoretical approach used to calculate the 
DR rate coefficient as function of electron collision energy 
was previously described in detail in Refs. [HI [131 US 112] ■ 
An intermediate step in the calculation is the following 
sum: 

^^,T.(^N + 1)P^V, (1) 

N 

N r 

where P^, ' is the probability of dissociative recombi- 
nation for a single e^+H;^ colhsional event. Here the 
system of the ion and the incident electron are described 
by the initial channel o' which describes both the initial 
rovibrational state of the ion |o') and its angular momen- 
tum coupling with the incident electron, assumed to be 
in a p-wave throughout. (See for instance, Eq. (4) and 



the corresponding discussion of Ref. [S].) The label F 
refers to the symmetry of the initial state of the system 
(the ion plus the electron at infinity), while and N 
are the angular momenta of the ion and neutral Hj]" + 
system, respectively. The quantum numbers T and N 
are conserved during the collision, along with the total 
energy E. 

The key ingredient in the calculation of the DR prob- 
N r • 

ability P^, ' is the energy-dependent scattering matrix 
S^^P{E), representing the scattering of the electron from 
the initial channel |o') to a final channel |o). The scat- 
tering matrix S'p'^^^(£') depends on the total energy E 
of the system, with the electron kinetic energy in chan- 
nel o' equal to Co' = E — Eq' where Eq' is the energy of 
the target ion in channel o'. In our theoretical treatment 
[81 I2TH23] . the complete energy-normalized wave function 
of the -I- e~ system at large distances r between the 
ion and the incident electron is written in the following 
way (see also Eq. (2.55) of Ref. [23) 



c 

exp{-ikor - i(l)o)Soo' - exp{ikor + i(j)o)S^^r{E) 



(2) 



The second sum above represents the scattering part of 
the wave function: The incident electron arrives in chan- 
nel |o') and can escape into one of the channels open for 
ionization. Those energetically-open ionization channels 
are labeled with an index o. In the above expression, 
we have explicitly written the electronic wave function 
of the incident electron only (r-dependence). The two 
other electrons of are considered to be frozen, but 
accounted for implicitly, because the scattering matrix 
(or quantum defect function used to calculate S^^'^^(E)) 
is obtained with those electrons present in the system, 
of course. The factor $0(0;) represents the ionic rovibra- 
tional wave function with lo the spatial and spin coor- 
dinates of the ion as well as the angular and spin coor- 
dinates of the scattering/Rydberg electron, i.e., includ- 
ing the angular and spin wavefunctions of the outermost 
electron and corresponding angular momentum coupling 
factors. The first sum in Eq. Q describes the part of 
the wave function in channels closed for ionization. The 
factor $c(w) represents the corresponding rovibrational 
wave function of the ionic core, again with the angular 
and spin degrees of freedom of the outermost electron. 
Wc^r^Vc) is the exponentially decaying electronic radial 
function in the closed channel, unity normalized, and A 
denotes an antisymmetrization operator for the electron. 
The coefficient Zco'{E) gives the relative contribution 
(amplitude) of the given closed channel to the total wave 



function: The relative probability to find the system in 
closed channel |c) is \Zco'{E)\'^, and this can be used to 
classify resonances. In a standard scattering theory for- 
mulation, an additional term would also be present to 
describe the dissociation channels of the neutral molecu- 
lar system. In our treatment, use of a complex absorbing 
potential in the vibrational degrees of freedom permits us 
to describe the loss of dissociative flux without explicitly 
computing individual scattering matrix elements into the 
dissociation channels. Lastly, in Eq. the phase factor 
(j)o denotes the usual Coulomb phase contributions, i.e., 
</>o = -lo 7r/2 + ao + (l/ko) ln2kor. 

Like the energy-dependent scattering matrix Sp^^^{E), 
the coefficients Zco'{E) are calculated from the smooth 
scattering matrix Si'i |24| . approximated here as energy- 
independent over an energy range 1-2 eV. First, the ma- 
trix 

Z,,,{E) = e-"''= {S,, - e-2'-c)-i S,,, , (3) 

is calculated. Here, Vc is the x diagonal ma- 
trix with nonzero elements v^. = Xj yj2{Ec — E) (effec- 
tive quantum numbers and atomic units); Ec is the ion 
energy for the closed channel |c). The matrix Z{E) cor- 
responds to energy-normalized electronic wave functions 
in closed channels. Therefore, to obtain amplitudes Z{E) 
for unity-normalized electronic functions, which corre- 
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spond more closely to standard spectroscopy classifica- 
tions, we have to multiply each element of Z{E) with 

3/2 

Z,,,{E) = v~''/^Z,o'{E). (4) 

Eqs. ^ and Q give the relative probability \Zco'{E)\^, 
proportional to the fraction of time spent by the electron 
in the closed channel c when incident in channel o', at 
any energy E. Because every channel (closed or open) is 
characterized with a specific ionization threshold energy 
Ec or Eo, the effective quantum number Vd for the domi- 
nant closed channel \d) is also obtained immediately. For 
energies E far from a resonance, several closed channels 
may have competing contributions. At a resonance, one 
closed channel usually dominates clearly, although for a 
"complex resonance" multiple closed channels can be im- 
portant for DR ^25, . 

IV. RESULTS AND DISCUSSION 

A. Imaging Results 

The 2D kinetic energy distributions of IR22p and 
Ipenn-09 are plotted in Fig. [T^). These distributions 
are projections of the total energy released, -E^g,., to the 
three H atoms upon dissociation, thereby revealing in- 
formation on the amount of internal energy available in 
the parent ion. The temperatures of the stored ion beams 
are derived with forward Monte-Carlo simulations, which 
describe the DR experiment from dissociation up to de- 
tection and analyze the simulated events like real data. 
The Dalitz plot presented in Ref. [7 is used to include the 
linear geometrical bias of the dissociation dynamics. An 
effective electron target-length is included to account for 
the overlap and merging and de-merging regions between 
the electron and ion beams. Using Eq. (6) of Ref. [TU], 
in both the measured and simulated data, the distance 
s '--^ 12 m from the center of the electron target to the 
detector is applied to give the tranversal energy E-j:^^ in 
terms of the inter-particle distances R. As the rotational 
population and the state-dependent DR rate coefficients 
are unknown, a Boltzmann distribution is modeled as- 
suming a rotationally independent rate. The rotational 
temperatures are derived from a fit between the simu- 
lation and the measurement using the temperature as free 
parameter. A temperature of 380 K could be determined 
for IR22p (see also Table [l]). A re-evaluation of Ipenn-09 
gave a roughly ten times higher temperature of 3250 K, in 
good agreement with the initially reported ~ 3500 K [7] . 
As can be observed in Fig. [T^) , the simulations and mea- 
surements deviate slightly, which is probably due to the 
approximations made. The shape of the simulated dis- 
tribution, however, comes much closer to the measured 
data then before, when the linear geometrical bias from 
the Dalitz data was not included. 




FIG. 1: (a) The measured 2D distributions ol the kinetic en- 
ergy released in the DR ol into H -I- H -1- H. In blue, 
the distribution for the beam produced with para-H2 in 
the 22-pole trap, IR22p, and in red, the distribution for the 
Hj}" beam produced with the Penning source, Ipenn-09. Also 
shown are the simulated distributions at 0, 380, and 3250 K. 
All distributions are normalized to unit area, (b) The semilog- 
arithmic zoom of the tail of IR22p (blue) and the previously 
measured I22p-09 (green). The simulations shown are at tem- 
peratures of 0, 200, and 400 K. 



The relative uncertainty, estimated from a 90% confi- 
dence interval of the fit, is about 35 K. Typically, the 
minimized was around 3. The systematic uncertainty 
is larger with ~ 135 K, dominantly determined by the un- 
certainty in the target length. In fact, the target length 
is estimated to be around 1.10 m, which is closer to the 
minimum than to the maximum overlap possible, giving 
rise to an asymmetric error bar (~iti3oK). The system- 
atic errors coming from the uncertainties of the ion-beam 
energy (~ 30 K) and kinetic energy release (~ 30 K for 
7 meV difference) have been significantly reduced as de- 
scribed below. 

The improved value for the H;^ dissociation energy 
[26H28] has led to an improved determination of the 
kinetic energy release in the H;^ DR. This kinetic en- 
ergy release has been recalculated to be 4.790 eV (from 
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iV+ = 1,K^ = 1) using the theoretical Hj^ dissociation 
energy that over the years has converged to 4.330 eV 
[53]. The experimental Hg dissociation energy of 4.337 
eV HZlllH], previously reported to be 4.373(21) eV P], 
gives only a 7 meV difference with theory, which yields 
about 30 K difference in the ion temperature. Addition- 
ally, the ion-beam energy for the measurements reported 
in Ref. [7 has now been analyzed in more detail, giving 
a more precise value of 4.096 ± 0.002 MeV. 

Figure [T|d) shows the tail of the 2D kinetic energy dis- 
tribution IR22p together with the previously measured 
imaging data of the DR of rotationally cool H|^, I22p- 
09. As can be seen, the distribution I22p-09 is colder 
than IR22p and even approaches the displayed 0-K sim- 
ulation between 4.8 and 5 eV, whereas IR22p follows the 
200-K curve more closely in this region. Above ^5.1 
eV, where the amplitude has dropped more than an or- 
der of magnitude already, the distributions approach the 
hotter curves, indicative of a non-Boltzmann character. 
Using the improved simulation and the updated energies 
to determine the temperatures of I22p-09 (shown) and 
I22n-09 (not shown), the minimum reduced with 90% 
confidence interval gives 200 ± 75K and 150 ± 50 K, re- 
spectively. Together with the systematic uncertainties of 
~ 135 K, the 22-pole ion-source temperature of 15 K 
cannot be entirely excluded. The re-evaluated tempera- 
tures are in agreement with the previous analysis stating 
<300K. The use of para-H2 or normal-H2 to produce the 
H;^ has no significant effect on the ion temperature of the 
stored beam. 

The main difference between IR22p and I22p-09 is the 
use of the electron cooler. The temperatures derived for 
IR22p (~ 380 K) and I22p-09 (< 300 K) are associated 
to the electron cooler on cooling energy for the entire 
storage time and 0.5 s only, respectively. Both measure- 
ments also differ in trapping times and gas pressures (see 
Table |l]) , however, these are not expected to yield dif- 
ferent temperatures. The lower (higher) He pressure is 
compensated for by the longer (shorter) trapping time 
and the trapping times should attain sufficient cooling 
power at the used buffer-gas densities as reported in Ref. 
[B| . This points to the electron cooler as source of heat- 
ing for the ion beam of the 22-pole trap. In Ref. [TU], 
at higher temperatures, the opposite effect is observed, 
where the electron cooler cools the ion-beam temperature 
of the supersonic ion source from 950 K to 450 K (see Ta- 
ble |l]). The heating of the 22-pole beam as opposed to the 
cooling of the supersonic beam by the electron cooler is 
suggestive of the electron cooler reaching an equilibrium 
between 380 and 450 K. The tranverse energy spread of 
the electrons in the cooler of ^ 12 nieV might be a factor 
in preventing cooling to lower temperatures. The tem- 
perature of the ion beam produced by the supersonic ion 
source is best approximated with a temperature of 950 K 
when the electron cooler is not used. Mutatis mutandis, 
the H;^ ion temperature after extraction from the 22- 
pole trap adheres to the < 300 K temperature and it is 
inconclusive whether any heating already occured in this 




FIG. 2: The H3 DR rate coefficients, R22n and Rpenn, as 
measured for the ion beam produced with the 22-pole trap 
and the Penning source, respectively, on a double-logarithmic 
scale. 



case. All rate coefficients, however, have been measured 
with the electron cooler turned on for the entire storage 
time. The rate coefficients measured with the 22-pole 
trap will therefore correspond to the heated temperature 
of ^ 380 K. Similarly, the rate coefficient measured with 
the supersonic ion source relates to the cooled tempera- 
ture of - 450 K. 

No imaging data exist for the rotationally hot H;^ beam 
from the Penning source with the electron cooler used for 
the entire storage time. Assuming the electron cooler is 
cooling at this high temperature, an ion-beam tempera- 
ture below 3250 K is believed to correspond to the DR 
rate coefficient of the highly excited beam. 



B. Rate Coefficient Measurements 

Figure [2] shows the high-resolution rate coefficients 
measured with high statistics as functions of the collision 
energy on a double logarithmic scale, before the reduc- 
tion performed in the detailed comparisons below. Ex- 
cept for the energy region of 2-8 meV, the rotationally 
hot H;^ rate coefficient Rpenn is larger than or equal to 
the rotationally cool rate R22n. Distinct structures are 
observed for the colder beam with the structures in par- 
ticular below 20 meV being more clearly resolved than 
ever before. The rate coefficient for highly excited ions 
is the highest resolution measurement for excited H;^ to 
date. Although much smoother, this rate coefficient still 
shows a remnant of the cold rate coefficient structures. 

The rate coefficients of the current and previous mea- 
surements (see Tabic |ll]) with rotationally cool H;^ are 
plotted in Fig. [Sj Here, the reduced DR rate coefficients 
are shown, which are the product {a)^/Ed and serve to 
eliminate the threshold dependence of the rate co- 

efficient a at low collision energies. All rate coefficients. 
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FIG. 3: The reduced DR rate coefficients from the rotation- 
ally cool Hjl" beam R22n compared with previously reported 
rates on a semi-logarithmic scale. For clarity, error bars are 
only shown for every 5th point. 



except R22n-09, are normalized to their 10-eV peak (not 
shown) with the absolute scale taken from Ref. 30. or 
Ref. lOJ in the case of Rsonic. R22n-09 has no 10-eV 
peak to normalize to. Previously, it was normalized to 
the 12-meV peak under the assumption that this reso- 
nance is temperature independent. As can be observed 
in Fig. |3]this may not be the case and the R22n-09 rate 
coefficient is instead normalized to the minimum around 
10 meV of R22n. Assuming the 10-meV normalization 
to be correct, the rate coefficients R22n and R22n-09 
look similar up to at least 20 meV. Above, comparison 
is difficult due to the large statistical errors of R22n-09. 
Both rates are expected to be similar as they have been 
measured under similar conditions. However, an even 
older measurement of the rate coefficient using the 22- 
pole trap as injector, R22n-05, exhibits less pronounced 
resonances around 6 and 12 meV, similar to the rate coef- 
ficient Rsonic from the supersonic ion source. The cause 
for the lack of structure in this region is presently not un- 
derstood. It is conceivable that the lower amplitudes of 
the 6- and 12-meV resonances is the result of hotter tem- 
peratures, similar to the resonance behaviour observed 
in the theoretical calculations discussed below. The ion 
temperature of Rsonic was estimated to lie around 450 K 
under the conditions of a rate measurement. This is 
slightly higher than the 380 K estimated to match R22n 
and R22n-09. 



C. Rate Coefficient Predictions 

Figure |4] displays theoretical rate-coefficients adapted 
from Ref. [8 for temperatures of 15, 100, 300, and 3000 K 
convolved with energy resolutions of AE„ = 25 /ieV and 
AE± = 500 fJ.eY, similar to experiment. The experimen- 
tal values given in Sec. [ll] are somewhat more conser- 



10= r 




FIG. 4; The Hg DR rate coefficients as predicted for 15, 
100, 300, and 3000 K, with the line thickness increasing with 
temperature. The 15 K rate is calculated using a 1:1 ortho- 
to-para ratio. 



vative, however, we do not expect this minor difference 
to be of importance. The relative para-/ortho-ratio of 
H;^ ions (including the multiplicity factor) used are in 
thermal equilibrium, rip/no = exp(A/fcBr)/2, except for 
the 15 K rate coefficient that assumes equal para and 
ortho populations (see Sec. |ll]). Using equal or ther- 
malized populations for higher temperatures makes lit- 
tle difference. The predictions show a generally higher 
rate for the hot H;^ than for the cold Hj^ , but with some 
local differences, as observed in the experiments. The 
rate-coefficient increase is most pronounced between 20 
and 200 meV, where the rate increases monotonically 
with increasing temperature. Above 200 meV, the rate 
coefficient decreases monotonically with increasing tem- 
perature. Three resonances around 9, 20, and 60 meV 
emerge with decreasing temperature, in between which 
there is no monotonic behaviour. It is expected that as 
the ion temperature decreases, the thermal averaging in- 
creases the contribution from the lower angular momen- 
tum states, and any individual resonance seen at low tem- 
perature should appear more strongly than in a higher 
temperature spectrum. Furthermore, the change in the 
amplitude of the rate coefficient seems to decrease with 
increasing temperature. Whereas the step from 15 to 
300 K changes the rate considerably, a generally smaller 
effect is observed when increasing the temperature from 
300 to 3000 K. 



D. Study of the Resonant Structure 

Figures [5^) and |6^) adopt a linear-linear plot to com- 
pare the experimental reduced DR rate coefficient of rota- 
tionally cool Hj^, R22n, with theoretical reduced DR rate 
coefficients in greater detail, up to a collision (detuning) 
energy of 0.06 and 0.5 eV, respectively. The theoretical 
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FIG. 5: (a) The reduced Hg" DR rate coefficient versus colli- 
sion (detuning) energy, as measured using the beam from the 
22-pole-trap (R22n), is compared with the theoretical rate 
coefficient calculated for a source of thermal ions at 300K, 
over the range from to 0.06 eV. The relevant excited rota- 
tional thresholds that can produce an infinite Rydberg series 
of np-resonances just below each of them are shown as ver- 
tical lines. The numerical values of those threshold energies 
and their corresponding quantum numbers are shown in Ta- 
ble |III[ For clarity, error bars are only shown for every 5th 
point, (b) Over the same energy range, theoretical spectra 
for the reduced DR rate coefficient are shown at four difi'er- 
ent temperatures (lOOK, 300K, 500K, 700K), with the line 
thickness increasing with temperature. 



curve is calculated using an assumed Hg ion temperature 
of 300 K, which is not far from the estimated ion beam 
internal temperature as was discussed above. Parts (b) 
of these two figures display the theoretically expected 
temperature-dependent rate coefficients from 100 K to 
TOOK in the corresponding low and high energy ranges. 

The following observations follow from studying Figs. 
[5]and[6l 

(1) Theory and experiment in the lower energy range 
from 0-0.06 eV exhibit clear disagreement in the shape 



FIG. 6: (a) The reduced Hj^ DR rate coefficient versus colli- 
sion (detuning) energy, as measured using the beam from the 
22-pole-trap (R22n), is compared with the theoretical rate co- 
efficient calculated for a sample of thermal ions at 300K, over 
the range from 0.05 eV to 0.5 eV. For clarity, error bars are 
only shown for every 5th point, (b) Over the same energy 
range, theoretical convolved spectra for the reduced DR rate 
are shown at four different temperatures (lOOK, 300K, 500K, 
700K), with the line thickness increasing with temperature. 



of the reduced DR rate versus energy, as is apparent in 
Fig. [5] The largest experimental resonance (A) appears 
to have no theoretical counterpart, and several promi- 
nent theoretical resonances (u, v, w, x, y, and z) have no 
discernable experimental companions. The lone region 
where theory resembles experiment in this region is the 
range from the local minimum (B) to a maximum (C), 
followed by a minimum (D). The theoretical local max- 
imum (c) is located about 1 meV higher in energy than 
the experimental maximum (C); it is possible that this 
approximate agreement is entirely fortuitous, given that 
this is the only feature in the lower energy range por- 
tion of the graph where the theoretical curve resembles 
experiment. 

(2) A working hypothesis concerning the low-energy 
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TABLE III: Ionization thresholds relevant to a p-wave elec- 
tron colliding at low energy with the Hj^ ion in difTerent initial 
ionic energy eigenstates characterized by either ortho (A2" or 
A2') or para (E" or E') symmetry. 



Target State {N+,K+) 


Upper State 


Ethreshold (meV) 


E" (1,1) 


(2,1) 


21.5 


E" (1,1) 


(3,1) 


53.4 


E" (2,1) 


(3,1) 


31.9 


E" (3,1) 


(4,1) 


42.0 


E" (4,1) 


(5,1) 


51.7 


E' (2,2) 
E' (3,2) 
E' (4,2) 
E' (4,4) 
A2' (1,0) 


(3,2) 


32.1 


(4,2) 


42.2 


(5,2) 


51.9 


(5,4) 


52.9 


(3,0) 


53.3 


A2" (3,3) 


(4,3) 


42.6 


A2" (4,3) 


(5,3) 


52.3 



discrepancy between theory and experiment in Figj5p) 
is that something problematic must be occuring in the 
treatment or the detection of the very highest Rydberg 
states in theory or experiment, respectively. It is in- 
triguing, and arguably more than a coincidence, that the 
4 rightmost theoretical resonances (centered near 0.02, 
0.03, 0.04, and 0.05 eV), none of which has any ex- 
perimental counterpart, are lying on and slightly below 
a rotational threshold supporting an infinite number of 
high Rydberg states. It should be kept in mind that 
the thresholds shown as vertical lines are total electron 
energies, whereas the theoretical and experimental DR 
curves are convolved with the asymmetric electron dis- 
tribution in velocity space, implying that the peaks in 
those curves correspond to higher total energies by ap- 
proximately 0.5 meV. If we turn the width of these reso- 
nances into a range of principal quantum numbers, these 
regions of enhanced DR are predominantly contributed 
by Rydberg states having high principal quantum num- 
bers, n > 50. A slightly suspicious point concerning the 
theory is that normally one would expect that a simple 
one-channel Rydberg series of p-states converging to an 
excited rotational threshold of the ion should decay by 
autoionization rather than by dissociating [31j. A high 
DR rate in this case is expected to occur only if there is 
a lower-n perturber, e.g. attached to a (01)^ state of vi- 
brational angular momentum, straddling that threshold 
energy range as well. On the other hand, we cannot rule 
out that theory might be correct in expecting enhance- 
ments at those rotational thresholds; in this case what 
should be explored is the possibility that the Rydberg 
states are either destroyed or their angular momentum 
quantum numbers are changed by external fields in the 
storage ring. 

(3) The higher energy range of Fig. [6] shows a good 
general agreement in magnitude and overall shape be- 
tween the experimental and theoretical curves, consis- 
tent with previous studies. The temperature dependence 
overall tends to casue a gradual increase with increasing 
internal temperature of the target ions, but in some lo- 



calized energy ranges exceptions are evident where the 
lower temperature gives a higher DR rate. 

(4) The DR rate coefficient in dependence on electron 
energy falls off approximately as l/\/Ed as the energy 
increases, over the range 0-0.3 eV, and at 0.3 eV (the 
first vibrational threshold as expected from the theoret- 
ical DR mechanism [Ml [23j [25]) the rate drops quickly 
from regions J and j to K and k, although somewhat 
more steeply in the theoretical treatments than in ex- 
periment. Both theory and experiment show a plateau 
region in reasonable agreement midway between features 
J and K. 

(5) The region of collision energies ranging from 
0.05 eV to 0.5 eV exhibits one of the best overall 
agreements between theory and experiment in the DR 
of H;^ . The shapes of both the experimental curve (from 
G-H-I-J-K) and the corresponding theoretical curve 
(ranging from g-h-i-j-k) show a level of similarity not 
observed before. In localized regions, however, theory 
is lower than experiment by 30-50%, especially near 
features G, H, and I. Resonant features I and i appear to 
coincide closely in position, although this feature does 
not involve just a single isolated resonance. Instead, 
the theoretical calculations in this energy range indicate 
many resonances from a number of different symmetry 
classes contributing to this structure. 

Figures [7^) and b) show a linear-linear plot of the ex- 
perimental reduced DR rate coefficients of the rotation- 
ally hot H;^, Rpenn, in comparison with theoretical re- 
duced DR rate coefficients at 1000 and 3000 K up to a 
collision (detiming) energy of 0.06 eV and 0.5 eV, respec- 
tively. The H;^ ion temperature of 3000 K is close to the 
derived ion beam internal temperature from the imaging 
experiments. As the temperature matching the rate co- 
efficient Rpenn is expected to be cooler due to the use 
of the electron cooler, the rate at an ion temperature of 
1000 K is shown as well. The following observations may 
be noted: 

(6) There is some similarity to be observed in the 
shape of the rate coefficient for the higher energy range as 
indicated by features F' to F and f ' to i'. All of these fea- 
tures are however a combination of resonances from dif- 
ferent symmetries and the mere occurence of resonances 
together with the differing amplitudes makes any defi- 
nite statement difficult. The experimental and theoret- 
ical rates are in good agreement in the energy region 
above J'(j')j '^i^^ the theoretical rate dropping slightly 
more steeply at the vibrational threshold around 0.3 eV, 
similar to observation (4) for the rotationally cool H;^. 

(7) Structures A to E are observed both for rotation- 
ally cool (Fig. [5| and hot H;[(Fig. [7)3) with the am- 
plitude of feature A increasing most with lower temper- 
ature. This could support the argument that the more 
pronounced structure observed in R22n and R22n-09 (see 
Fig. [3]) is due to a lower rotational ion temperature. 
There are several possible explanations for the consis- 
tent shape of the A-E structure at varying rotational 
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FIG. 7: (a) The reduced Hj^ rate coefficient versus coUision 
(detuning) energy, as measured using the beam from the Pen- 
ning source (Rpenn), is compared with the theoretical rate 
calculated for a source of thermal ions at 1000 K (purple) and 
3000 K (red), over the range from eV to 0.06 eV. For clarity, 
error bars are only shown for every 5th point, (b) The same 
for the energy range of 0.05 eV to 0.5 eV. The line thickness 
increases with temperature. 



temperature. One explanation might be that the A-E 
structure is a consequence of higher populated low iV+ 
levels than assumed for a lO'^ K Boltzmann distribution. 
The rotational population of the beam is probably 
not Boltzmann-distributed. Additionally, higher rate co- 
efficients for low 7V+ levels than for high A''+ levels would 
lead to a higher contribution from the lower levels. This 
could give rise to an apparent remnant cold structure in 
the hot H;^ rate coefficient, as can also be observed in 
Fig. [2j Another explanation is that the seemingly in- 
variable behaviour is fortuitous. Features A-E consist of 
multiple resonances and are most probably contributed 
by low Rydberg state perturbers (n < 9). The contribu- 
tions of low and high levels will therefore give rise to 
resonances at similar energies, producing similar though 
different features. This is well demonstrated in regions 



t' and u' where the shapes of the resonance features for 
1000 and 3000 K are similar, but where in this case the 
different origins of the features are revealed by a clearly 
visible change in position. 

(8) A larger discrepancy appears in the overall ampli- 
tude of the rate coefficient than for the low-temperature 
counterpart. The average value of the predicted reduced 
DR rate is clearly lower than the one observed, in both 
energy regions. Enhancement of the rotational temper- 
ature does not narrow this gap as it did in the low- 
temperature case. Additonal rotational states have little 
more effect on the total rate as can be observed by the 
similar average rates for 1000 and 3000 K temperatures. 
It should be kept in mind, however, that the 3000 K theo- 
retical rate coefficient is not fully converged with respect 
to the highest ionic angular momenta in the Boltzmann 
distribution. 



V. CONCLUSIONS 

The temperature derived by the imaging measure- 
ments for the 22-pole trap under the conditions of a rate 
measurement is clearly higher than had originally been 
expected. As suggested in the theoretical work of Ref. 
8J, a higher temperature of the stored ions has to be 
assumed for the measured rate coefficients of the rota- 
tionally cool H^. Although Ref. [8_ proposed an ion 
temperature of 1000 K for the storage-ring experiments, 
it seems that 300 K is sufficient to attain a similar aver- 
age rate. Disagreements up to an order of magnitude are 
hereby reduced to a factor of roughly 2. 

The experimental work reported in Ref. [10] suggests 
that the extraction from the supersonic ion source could 
be the cause for the heating of the ion temperature (up to 
950 K) . It is unclear whether ion extraction causes heat- 
ing also in the 22-pole trap. The ion-beam temperature 
derived from previous imaging measurements, I22p-09 
and I22n-09 (see Table |l]), is confirmed to be < 300 K 
using an improved simulation of the dissociation events 
and including an updated kinetic energy release and ion 
beam energy. With the statistical and systematic uncer- 
tainties amounting to roughly 135 K, the trap tempera- 
ture of '-^ 15 K cannot be entirely excluded. Heating is, 
however, observed inside the ring. The < 300 K temper- 
ature applies to conditions where the electron cooler is 
used for only 0.5 s. When employing the electron cooler 
for the entire storage time, as was the case for IR22p, the 
ion temperature is raised to ~ 380 K, implying heating 
by the electron beam of the cooler. This is the opposite 
effect as was observed in Ref. [TU] for the supersonic ion 
source, where the electron cooler seems to cool the H;^ 
beam from initially ~ 950 K to ~ 450 K. The cooling of 
hot beams and heating of cold beams may be an indica- 
tion of an equilibrium being attained around 400 K. If 
so, the rotationally hot H;^ beam of the Penning source 
will be colder than the 3250 K temperature derived by 
imaging (Ipenn-09) with the electron cooler on for only 
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0.5 s. 

The rate coefficients measured at the storage ring 
TSR are all obtained with the electron cooler on cooling 
conditions for the entire storage time, implying ^ 380 K 
with the 22-pole trap and ^ 450 K with the supersonic 
ion source [10]. Presently no rate-coefficient measure- 
ment with a confirmed temperature below 300 K exists. 
The present work suggests that in order to obtain a rate 
coefficient measurement at lower ion temperatures, the 
electron beam of the electron cooler should be excluded 
from the measurement or its influence reduced by mini- 
mizing its transverse energy spread. The electron target 
having an electron beam with a much smaller velocity 
spread (1 meV) could be used for phase-space cooling 
the ion beam instead. Its toroidal regions are also de- 
signed to (de)merge the electron beam faster, reducing 
any possible temperature change from non-zero energy 
collisions. To eliminate possible heating from ion-source 
extraction, which so far can be neither confirmed nor ex- 
cluded, the helium buffer-gas may be pulsed into the trap 
during initial cooling only, reducing the particle density 
upon extraction. 

At the currently estimated temperature of ^ 380 K, 
95% of the ions are distributed over 15 rotational states 
up to — 5 instead of the previously reported cold 
rovibrational populations of predominantly = 1. 

These are too many rotational states to make a detailed 
study of the contributions of specific symmetries. DR 
measurements of real rotationally cold H;^ ion beams are 
necessary to advance further. The level of agreement 
between theory and experiment is however encouraging, 
especially in the energy range 0.06 eV to 0.5 eV, where 
a number of features that modulate the overall energy- 
dependent DR rate coefficient appear to be mirrored in 
theory and experiment. This comparison is the most de- 
tailed one carried out to date. Nearly all of the major 
discrepancies at collision energies below 0.06 eV arise in 
a narrow range approximately 4 meV below ionization 
thresholds supporting very high Rydberg states. This 



points to the need of combined efforts by theory and ex- 
periment in the range of very high Rydberg principal 
quantum numbers n > 50 in order to track down the 
origin of this major disagreement. 

The present analysis stresses the need for continuing 
efforts to develop an in-situ method for determining and 
if possible controlling the rotational state distribution of 
the ions in the storage ring. Preliminary work on an 
in-situ probing procedure is being performed |32j . Addi- 
tionally, the predicted 10:1 higher rate coefficient for the 
para- (1,1) than for the ortho-(l,0) state can be studied in 
more detail when low temperatures are attained. Further 
improvements to the theory which are under considera- 
tion elsewhere are (i) the elimination of the currently- 
implemented rigid-rotator approximation for the ro- 
tational degrees of freedom; (ii) inclusion of I — mixing 
effects due to the non-spherically-symmetric interaction 
experienced by an electron that collides with the ion; 
and (iii) continued exploration of external field effects on 
Rydberg-mediated DR. 
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